Sorption of proteins to charged microgels: characterizing binding isotherms and 

driving forces 
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We present a set of Langmuir binding models in which electrostatic cooperativity effects to protein 
sorption is incorporated in the spirit of Guoy-Chapman-Stern models, where the global substrate 
(microgel) charge state is modified by bound reactants (charged proteins). Application of this 
approach to lysozyme sorption to oppositely charged core-shell microgels allows us to extract the 
intrinsic, binding affinity of the protein to the gel, which is salt-concentration independent and 
mostly hydrophobic in nature. The total binding affinity is found to be mainly electrostatic in 
nature, changes many orders of magnitude during the sorption process, and is significantly influenced 
by osmotic deswelling effects. The intrinsic binding affinity is determined to be about 7 ksT for 
our system. We additionally show that Langmuir binding models and those based on excluded- 
volume interactions are formally equivalent for low to moderate protein packing, if the nature of the 
bound state is consistently defined. Having appreciated this, a more quantitative interpretation of 
binding isotherms in terms of separate physical interactions is possible in future for a wide variety 
of experimental approaches. 



I. INTRODUCTION 

Functionalized colloids and nanoparticles play an in- 
creasingly prominent role in the development of bioma- 
terials and imminent biotechnological applications, for 
example, drug delivery, enzyme biocatalysis, or control 
of gene expression^— In particular, multiresponsive hy- 
drogels are of great interest due their biocompatibility, 
resemblance to biological tissue, and tunable viscoelas- 



tic properties 
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Dispersed in water these colloidal 



microgels create an enormous surface that may be taken 
as a model for soft biological interfaces. Hence, protein 
storage, activity, and uptake properties may be changed 
by physiological stimuli, such as pH, salt concentration, 
and temperature. Moreover, colloidal hydrogels qualify 
for a number of applications, e.g., as drug carrier de- 
vicesi^r— However, the detailed control of functionality 
requires a quantitative understanding of the underlying 
physical interactions between microgels and biomolecules 
in the aqueous environment. 

Recent studies have demonstrated that protein sorp- 
tion to nanoparticles is mostly driven by global, nonspe- 
cific electrostatic interactions and more local, probably 
hydrophobic interactions i^ 7 - ! 12 ' 19- — The balance between 
those two is highly system-specific and can be manipu- 
lated by chemical functionalization or copolymerization. 
Charged microgels, for instance, can be used to favor or 
disfavor the sorption of net-charged proteins, while their 
osmotic swelling and storage volume can be tuned by 
pH, salt and charge densit y 10 ' 11 ^ 5 ' 27 essentially via the 
Donnan equilibrium. 2 - However, during the uptake of the 
charged, polyionic proteins, swelling and Donnan equi- 
libria are typically changing in an interconnected fash- 
ioni 19 ' 24 " — These highly cooperative effects render the 
interpretation of binding isotherms, and thus the sep- 
aration and quantification of global electrostatic and lo- 
cal hydrophobic contributions to binding, a difficult task. 



Moreover, binding affinities in these systems depend on 
protein load which presents an additional complication 
when modeling the adsorption isotherm. 

Up to now, the modeling of protein uptake has 
been done often by using the standard Langmuir 
isotherm. ) 29 i 30 in particular when evaluating protein ad- 
sorption as measured by isothermal titration calorimetry 
(ITC^i^i^i^MiL In the standard Langmuir approach, 
however, protein association with single, independent 
binding sites is assumed which neglects electrostatic co- 
operativity effects and volume changes during sorption. 
Additionally the term 'binding' of proteins to soft poly- 
meric layers and hydrogels is somewhat ill defined, as the 
system may remain in a fluid-like state where proteins are 
still mobile on average, albeit slower than in bulkj 13 i 32 
Consequently the stochiometry and binding affinities to 
'sites' in the gel obtained from Langmuir fitting are not 
so easy to interpret. 

Alternatively, the hydrogel matrix may be viewed as a 
homogeneously charged background to the mobile pro- 
teins where saturation of sorption may set in due to 
excluded- volume (EV) packing. In recent binding models 
based on such a viewj 25 i 33 ~ — it is typically assumed that 
proteins can be treated as simple charged hard spheres, 
i.e., hard polyions. The electrostatic problem can then be 
tackled by approximative Poisson-Boltzman (PB) 'cell' or 
'box' models, — ~— where electrostatic cooperativity and 
osmotic ion effects to deswelling can be included. So far 
it has not been attempted, however, to separately treat 
more specific effects in EV-based binding models such 
as hydrophobic interactions or restraints of the configu- 
rational protein degrees of freedom in the bound state. 
Also the relation between Langmuir and EV models, if 
any, is unclear. However, modeling protein adsorption by 
a meaningful and physically sound isotherm is the pre- 
requisite for a quantitative understanding of the driving 
forces. 
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Here we present an in-depth discussion of the isotherms 
suitable for modeling protein sorption into microgels and 
soft polymeric layers in general. This discussion will 
then provide the basis for a detailed investigation of 
the driving forces of protein binding. We demonstrate 
that electrostatic cooperativity and the effects of mi- 
crogel volume changes can be introduced into standard 
Langmuir models following the spirit of Guoy-Chapman- 
Stern theory for binding of charged molecules to charged 
surfaces ) 46 We test the performance of the binding- 
models by fitting to binding isotherms obtained previ- 
ously from ITC of chicken egg white lysozyme sorption 
onto a negatively charged coreshell microgel* 2 - 3 - The core- 
shell particles consist of a polystyrene core onto which 
a charged poly (N-isopropylacrylamide-co-acrylic acid) 
(NiPAm) network is attached. Previous characterization 
of this system shows that it is an ideal model system as it 
reaches full equilibrium with high binding affinities, and 
the globular lysozyme maintain its folded state in the 
microgel with even enhanced activity^ 3 - We then demon- 
strate that we can separate out the global electrostatic 
contribution leading to consistent values for the salt con- 
centration independent binding affinity Finally we show 
that EV-based binding models are formally equivalent to 
the Langmuir approach in the low-packing regime if the 
bound state is consistently defined. Consequences to the 
interpretation of Langmuir models are discussed. 



II. EXPERIMENTS: MATERIALS AND 
METHODS 

A. Materials 

In this study the same batch of microgel dispersion was 
used as in previous workj^ In brief, the polystyrene core 
was synthesized by emulsion polymerization in the first 
step. After purification of the core particles, the microgel 
shell, containing 5 mol-% N,N-methylenebisacrylamide 
(BIS) crosslinkers and 10 mol-% acrylic acid with re- 
spect to the amount of NiPAm, was polymerized on the 
polystyrene core by seed polymerization. After purifica- 
tion the particles were transferred into buffer solution by 
ultrafiltration against 10 mM MOPS buffer at pH 7.2. 
In this preparation state the gel is swollen and strongly 
hydrated with more than 90% volume fraction of water. 
Dynamic Light Scattering (DLS, Malvern Instruments) 
was applied to determine the hydrodynamic radius of the 
polystyrene core to i? CO ro = 62.2 ± 0.7 nm and the radius 
of the total core-shell microgel to R ~ 129 to 172 nm 
depending on the solution condition. 23 

The ITC experiments were performed using a VP-ITC 
instrument (Microcal) as described previously >2S Briefly, 
a total of 300 ml of lysozyme solution (0.695 mM) was 
titrated into the sample cell filled with 1.4 ml of buffer- 
matched microgel dispersion at a c TO = 8.42 • 10~ 7 mM 
concentration. The pH-value was held constant at 
pH=7.2. The experiments were performed at 298 and 



303 K. The following three buffer systems were used: (i) 
10 mM MOPS, 2 mM NaN 3 (7 mM ionic strength); (ii) as 
(i) with additional 10 mM NaCl (17 mM ionic strength); 
(hi) as (i) with additional 25 mM NaCl (32 mM ionic 
strength). The incremental heat changes AQ were mea- 
sured during the course of the titration experiment where 
the protein solution is stepwise injected into the microgel 
dispersion. The integrated heat change after each injec- 
tion was corrected by the heat of dilution. 

III. THEORY 
A. Basic model 
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FIG. 1: Illustrative sketch of the model for one core-shell mi- 
crogel in a protein dispersion. The microgel is represented by 
a sphere with radius ii between 129 and 172 nm including a 
core with radius -R CO rc = 62.2 nm. The gel with volume V has 
a mean charge density c a = N g /V . Protein and monovalent 
salt concentrations in the bulk region are c p and c s , respec- 
tively. The proteins have radius R p and valency z v . Nb p is 
the number of bound proteins. 

In our model one core-shell hydrogel particle is mod- 
eled as a perfect sphere with radius R having a core with 
radius i? C ore as depicted in Fig. 1. The net gel volume is 
therefore V — 4ir(R 3 — i? 3 oro )/3.0. Potentiometric mea- 
surements show that a total number of N g = 4.9 • 10 5 
charged monomers are present in the gel. The network 
monomers are assumed to be homogeneously distributed 
in the gel as justified by small angle scattering^ The 
mean gel charge density is z g c g = z g N g /V, where z g is 
the monomer charge valency. Since we work at a pH 
~ 7.2, much larger than the pK a -value of ~ 4.6 of poly- 
acrylic acid, and deal with weakly charged gels (with a 
charge fraction ~ 1/10), we can safely assume that charge 
regulation plays a minor role and z g = — 1 in the follow- 
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ing. The salt is monovalent with a bulk concentration 
c s = c + = c_ , where c+ and c_ are the concentrations of 
cations and anions, respectively. The Bjerrum length is 
Ib = e 2 /(47reeofcsT), where e is the elementary charge, 
eo and e(T) the vacuum and bulk permittivities, respec- 
tively, and ksT = the thermal energy. In our sys- 
tems Ib — 0.7 nm, while the electrostatic Debye-Hiickel 
screening length k^ 1 — (SttIbCs)^ 1 ^ 2 in bulk is between 
~ 3.7 nm (7 mM salt), ~ 2.4 nm (17 mM salt), and ~ 1.7 
nm (32 mM salt). 

The lysozyme proteins are modeled by monopolar, ho- 
mogeneously charged spheres with diameter a p = 2R p , 
assumed to have a valency z p = +7 (i.e., a net charge of 
+7e) estimated from titration^ 8 - Hence, multipolar con- 
tributions, e.g., dipolar or quadropolar or charged-patch 
interactions to sorption are neglected. The proteins are 
found in bulk solution with density c p . No dispersion 
or other attractive nonelectrostatic interactions between 
the proteins are considered, and we assume a homoge- 
neous distribution of proteins in the gel region^ Aggre- 
gation of proteins in the gel is unlikely because at full 
load the system is still below the solubility threshold of 
a bulk system at comparable pH and electrolyte concen- 
tration^ The number of bound protein inside the gel is 
denoted by Nb p giving rise to an internal protein pack- 
ing fraction r\ = Nb p Tra p /6V. Finally, the gel and the 
aqueous buffer are modeled as a continuum background 
with elastic, dielectric, and osmotic properties as detailed 
in the following sections. The introduced quantities and 
variables are summarized in Tab.l. 



B. Electrostatic contributions to the interactions of 
the proteins with the charged gel 

1. Donnan equilibrium 

Homogenization of the charged gel by salt ions leads 
to a electrostatic potential difference between bulk and 
gel, the Donnan potential. The latter can be derived 
by assuming electroneutrality in the gel for monovalent, 
ideal ions leading to2£ 



variable 


meaning 


R 


gel radius 


V 


gel volume 


c s 


charged monomer concentration 


c s 


bulk salt concentration 


Kg 


inverse Debye length in gel 


K b 


inverse Debye length in bulk 


Cp 


bulk protein concentration 


tot 
t,p 


total protein concentration 


X 


molar ratio x = c P Qt /c m 


Nip 


number of proteins in the gel 


N 


number of Langmuir binding sites 


e 


protein load = N bp /N 




effective protein diameter 


Rp 


effective protein radius 


n 


protein packing fraction in gel 




Donnan potential 


K 


binding constant 


K m 


gel modulus 


constants 


meaning 


Rcorc — 62.2 nm 


core radius 


N g = 4.9 ■ 10 5 


number of charged monomers 


Zg = -1 


monomer charge valency 


Zp — 7 


protein charge valency 


c m = 8.42 ■ KT 10 M 


microgel concentration 


Ib m 0.7 nm 


Bjerrum length 


wo=l/mol 


standard volume 



TABLE I: Frequently used variables and constants. 



mM, thus \A<j>\ is on the order of unity (~ 25 mV) and sig- 
nificant for the interaction with charged macromolecules. 

The Donnan equilibrium leads to an osmotic pressure 
difference 7r; on between gel and bulk ions and, in the ideal 
gas limit, is given by the difference of internal and exter- 
nal ionic concentrations 28 



faion = C s e- A * + 



= A0 



2c, 



(2) 



2c s vV + 1 - 1 = 2c s [cosh A(f> - 1]. 



For the salts and concentrations considered in this work 
corrections due to nonideal activity are small^ and can 
be safely neglected. 



ee ePA(f> = In 



+ V? 



1 



(1) 



2. Electrostatic protein-gel interaction 

The transfer of a charged protein from the bulk so- 
lution with salt concentration c s into the charged gel 
with monomer charge concentration c g is accompanied 
by a (Gibbs) transfer free energy change per particle (or 
chemical potential). If we assume the gel monomers to 
be mobile, this can be modeled by the transfer of a spher- 
ical polyion with radius R p and valency z p from one bath 
at salt concentration c s and zero potential to another at 
mean potential A<f> and salt concentration c g . The dif- 
state. Given bulk salt concentrations between 7 and 32 ference in solvation free energies on a Debye-Hiickel level 



where y = z g c g / (2c s ) is the charge ratio between gel and 
bulk charge densities and the tilde in A</> denotes the 
dimensionless potential scaled by e/3. 

We estimate a mean separation of about 2.2 nm be- 
tween two charged monomers on the same polymer 
within the gel, considerably larger than the Bjerrum 
length I b — 0.7 nm. Thus, charge regulation effects by 
counterion (Manning) condensation or inhomogeneity ef- 
fects can be safely neglected^Sr— For our microgel we 
further estimate the charged monomer concentration to 
be c g = N g /V ~ 40— 100 mM, depending on the swelling 
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is then 



f3AG ei = z p A(t>(y) 



ZplB 

2R n 



+ KgRp 



K b Rp 
1 + K b R, 



(3) 



where we defined the inverse screening length n g — 
y^8nlBc g in the gel, characterizing approximately the 
screening by mobile polyelectrolyte charges and their 
neutralizing counterions (see Appendix A). The Donnan 
potential ([lj must now be corrected for the change of 
total net charge by the charge of bound protein via 



V = ( Z 9 C 9 + z pN bp /V)/(2c s ). 



(4) 



The first, purely entropic term in eq. ([3]) simply expresses 
the electrostatic transfer energy of a charge z p from bulk 
to a region at potential A(j)(y). The second term reflects 
the difference in the Born solvation free energies in a 
salty medium.— The first term is attractive if microgel 
and protein have opposite net charge, otherwise repul- 
sive, while the second one is attractive if c g > c s , other- 
wise repulsive. Since we deal with a weakly charged gel 
additional entropic effects from the release of condensed 
counterions can be neglected in our system^ 

In Appendix A a more rigorous derivation of eq. ([3]) 
is presented in the framework of a linearized PB cell 
model^— corrected for particle-particle interactions. 
While eq. ([3]) is valid strictly only in the limits of small 
salt concentrations and small protein charges, its quanti- 
tative validity seems to cover a wide range of protein load 
as long as the protein valency is not too high, as discussed 
in Appendix A. For our system with lysozyme (surface 
charge density ~0.17e/nm 2 ) in monovalent salt and the 
weakly charged gel (charge fraction 1/10), Debye-Hiickel 
level assumptions are justified. 



C. Osmotic and elastic swelling equilibrium 

In equilibrium the gel size is determined by a mechan- 
ical balance between osmotic and elastic pressures in the 
gel via^^i^-^ 



+ Kel = 0. 



(5) 



The osmotic term Tr osm has two contributions one from 
the solvent, the other from the ions. The first is often 
expressed by a de Gennes-like scaling type of relation 
oc V~ n , where typically n = 9/4 for neutral polymers 
under good solvent conditions . 58 ' 61 For charged networks 
corrections may arise^ while for our case of weakly 
charged and weakly stretched gels the 9/4 exponent likely 
to be valid i 57 ' 58 ' 62 (According to definitions by Dobrynin 
and Rubinstein,— we estimate an elecrostatic blob size 
of about 1.3 nm using a monomer length 0.35 nm and 
our charge fraction 1/10. This is indeed smaller than the 
correlation length of 3.2 nm given by our mean mesh size, 
estimated from the linker density.) 

The ionic osmotic pressure is dominated by the ideal 
gas pressure of ions in the gel as given by eq. ^ . This is 



counterbalanced by the elastic pressure n e i as represented 
essentially by the shear modulus G. Weakly charged net- 
works obey Gaussian chain statistic o 57 ' 59 , and G scales 
via the power law 7r e ; cx V~ m with m = 1/3. Various 
experiments corroborate that scaling] 57 ' 58 ! 62- — The ex- 
pression for the total pressure is 



TT = AV~ 

= AV- 



BV- m 



Kioniy) 



(6) 



where A and B (or Vq) are volume- independent con- 
stants. For high salt concentration the ionic contribu- 
tion vanishes, 7Tj „ = 0, and the equilibrium gel is in the 
'neutral' reference state with volume Vq. The bulk mod- 
ulus of the gel in this state is defined via /3K m (Vo) — 
-VdpP/dV\ T ,v = (23/12)A/Vo™. By fitting eq. © to a 
variety of salt concentrations (without proteins) we will 
obtain both the unknowns A and Vq and thus K m (Vo)- 

In first order the effect of protein addition to the gel is 
lowering of the gel net charge and the inhomogenization 
of the charge and electrostatic potential distribution in- 
side the gel. As shown in the PB cell model in Appendix 
A those effects can be approximately included in © by 
replacing n ion (y) with 

nf on (y) ~ w ion (y) + 2k B Tc s [4>{R c ) - A$\ sinh(A^), (7) 

with y = (zgCg + z p N bp /V)/(2c s ). The potential <p(R c ) 
is the electrostatic potential at the cell boundary R c in 
the PB cell model (Appendix A). 

Estimating other protein-induced contributions to the 
pressure from more local effects on elasticity, such 
as cross-linking by local binding^ conformational re- 
straints of the polymer network, or possible contributions 
from the protein osmotic pressure is challenging due to 
the lack of precise knowledge of the nature of the bound 
state and is out of scope of this paper. 



D. Binding isotherms 

The transfer of one protein into the gel will be accom- 
panied by a release of the binding enthalpy AH. If N bp 
proteins bind and we assume that the heat per protein 
does not change with load, the total heat released is 



Q(N bp ) = AHc m V tot N b p, 



(8) 



where Vtot is the total titration volume. In the ITC ex- 
periments Q is measured vs. the total protein concen- 
tration cp ot . Introducing the molar ratio x = cp ot /c m we 
can write Q(x) = AHc m Vt tNb p (x). The incremental 
heat Q'{x) = dQ/dx per molar concentration of protein 
is more sensitive to fitting and reads 



y(x)/(V tot cp ot ) = AHN' bp {x)/x. 



(9) 



For large binding constants and small x, almost all of the 
proteins immediately get sorbed, and N bp (x) ~ x, and a 
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plateau of heigh Q' (x) / {V t otC P ot ) ~ AH is expected, in- 
dependent of the binding constant. For large x, typically 
N bp (x) saturates and Q'[x) oc N' bp {x) = 0. Thus we rec- 
ognize that fitting to binding models is most sensitive 
to intermediate values of the molar ratio x, in the pre- 
saturation regime, near the inflection point of Q'(x). 



1. Standard Langmuir Binding Model 

The standard Langmuir model is based on identifying 
the association reaction A+B— >AB, with a binding con- 
stant K = [AB] / [A] [B] , where the square brackets denote 
concentrations. The basic assumptions in the Langmuir 
model are that ideal particles are binding to a fixed num- 
ber N of identical and independent bindings sites. The 
equilibrium constant in the Langmuir framework is usu- 
ally written a o 29 ' 30 



K 



e 



(l-e)cp' 



(10) 



where 9 = N bp /N denotes the fraction of bound protein 
Nt p and total sites N (see also derivation in Appendix B 
for the canonical ensemble) . The protein density outside 
the gel can be expressed by the total protein density in 
the sample minus bound protein density by c p = c* ot — 
NQc m - If it is assumed that N, K, and AH are protein 
concentration independent, solving ([ID| for 6 gives the 
total heat 



Q{x) = ^NAHc m V tot 



(11) 



with £ = 1 + x/N + \/{NKc m ). The fitting of 
Q' (x) I (V t otC P ot ) to the experimental data then yields the 
unknown constants K, N, and AH. Typically a sig- 
moidal curve is measured for Q'(x), where AH describes 
the plateau for the first injections (small a;), N the in- 
flection point, and K the sharpness of the transition at 
x = N. Fitting to Langmuir isotherms is thus most sen- 
sitive close to x ~ A, when the molar ratio equals the 
number of available binding sites. 

The Gibbs binding free energy in the Langmuir model 
is defined by 



pAG = -HK/vo), 



(12) 



where vq is the 'standard volume' which describes the 
'effective' configurational volume in one of the binding 
boxes (see also Appendix B). Thus the absolute value of 
AG depends on the magnitude of the standard volume vo, 
an often overlooked fact in literature.— Typically vo — 
1/mol ~ 1.6 nm 3 is chosen which is a reasonable choice 
for molecular binding where spatial fluctuations are on a 
nanometer length scale. While for quantitative estimates 
thus precise knowledge of the nature of the bound state 
is necessary, we will satisfy ourselves in this work with 
the standard choice vq = 1/mol. 



2. Extended Langmuir Binding: Separating 
Nonelectro -static and Electrostatic Contributions 

When many charged entities bind to charged regions a 
cooperativity effect comes into play due to the change of 
the global electrostatic properties during loading. This 
has been appreciated in the Guoy-Chapman-Stern the- 
ory for the binding of charged ligands to charged sur- 
faces4i~— Consequently, the total binding constant K = 
K(x) has to be defined more generally and split up into 
an 'intrinsic' part Kq and an electrostatic part via 



K{x) = K exp[-f3AG cl {x)} = 



Q(x) 



[l-Q{x)]c p 



(13) 



i.e., AG(x) — AG e i{x) + A Go- In our work we as- 
sume AG e \(x) to be given by eq. ([3]). The intrinsic, 
x-independent binding constant defines the intrinsic 
adsorption free energy 



pAG = -]n(K /v ), 



(14) 



which only contains contributions from specific interac- 
tion between the protein and the gel environment, and 
the nonspecific, leading order electrostatic effect has been 
separated out. Specific interaction may include local sol- 
vation effects - where the hydrophobic effect is naturally 
one of the biggest contributors - and possibly specific 
local interactions, such as salt bridges. 



3. 'Excluded Volume' (EV) Binding Isotherms 

An alternative binding isotherm that may be used for 
evaluating the ITC data is based on the equivalence of 
chemical potentials in bulk and inside the gel, where in 
addition to electrostatic and intrinsic binding effects the 
EV interaction between the (hard spherical) proteins in- 
side the gel are taken into account i 25 > 33 ~ — One assump- 
tion is here that in the very contrast to Langmuir mod- 
els, the proteins are not 'condensed' to fixed sites in the 
gel but are able to freely move around, under the re- 
straint only that their translational freedom is confined 
by packing (excluded volume). This ansatz yields the 
Boltzmann-like equation for the density of hard spheres 
in the gel 

N bp /V = C'c p exp(-/3AG ei - /3AG ) exp(-/3/i C s), (15) 

where £' is the partition function of the bound state more 
precisely defined later, and fics is the Carnahan-Starling 
(CS) excess chemical potential^ 



8r] - 9r, 2 + 3r, 3 

(1 ~ v) 3 



(16) 



representing the free energy of transferring one hard 
sphere to a solution with the packing fraction rj = 
(Ab p /l/)7T(7 3 /6. The effective diameter of the sphere 
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(j p = 2R p is expected to be close to the diameter of 
gyration or hydrodynamic diameter of the protein (i.e., 
Op — 3 — 4 nm for lysozyme)^ and will serve as a fit- 
ting parameter in the following. Since packing effects 
play a role in the EV approach we correct the gel volume 
by the polymer volume fraction. A PNIPAM monomer 
has an excluded volume roughly of 0.3 nm 3 . We es- 
timate 3. 7 TO 6 monomers in one microgel what makes 
Vpoiymor — 1.1-10 nm 3 . This corresponds to a polymer 
volume fraction in the range of 5% (pure, unloaded gel) 
to 12% (fully loaded gel in 7mM salt). 

In (|15l) we have also separately described the non- 
specific electrostatic part (AG e ;) and the intrinsic part 
(AGo) to binding as in the Langmuir model above. In 
the EV model, the contribution to AGo can be still 
thought of being induced by hydrophobic interactions or 
salt bridges but only in a weak-interaction sense, such 
that the particles in the gel are still translationally free. 
However, we additionally consider £' in eq. (1151) . the par- 
tition function of the protein in the bound state, which 
can include contributions from vibrational and orienta- 
tional restrictions of the protein's degrees of freedom 
within the gel, e.g., by partial sliding on polyelectrolyte 
chains iZS 



4- Equivalence Between the Langmuir and the EV 
Approach 

For small protein packing fractions we may linearize 
the CS chemical potential in eq. (fT5")) . We obtain 



N bp 
V 



C'c p exp(-/3AG ei - f3 AGo) (1 - 2B 2 N bp /V) (17) 



where we identified the second virial coefficient of hard 
spheres B 2 = 27rcr 3 /3. Physically 2B 2 describes the vol- 
ume excluded to the centers of the other spheres taken 
by one sphere. We rearrange to obtain 



exp(-/3AG ei - f3AG ) = 



N bp /V 



(1 - 2B 2 N bp /V)c p ' 



If we now make the substitution N = V/ {2B 2 ) and 9 = 
Nbp/N as in the Langmuir model above, we end up with 
the standard Langmuir form (flU)) 



C'2B 2 exp(-/3AG ei - /3AG ) = 



e 



(l-e)cp 



(19) 



If £' = 1, thus the EV treatment is equivalent to a Lang- 
muir picture where a bound ideal gas particle has a con- 
figurational freedom (volume) of 2B 2 , i.e., it can freely 
move around in one of the N binding boxes. However, 
as discussed above, in the Langmuir-type bound state 
the configurations are restricted to an effective configu- 
rational volume vq with respect to 2B 2 , such that the 
partition function £' = vo/(2B 2 ). With that definition 
we exactly end up with the standard Langmuir model 
eq. (29- 



Thus, in the approximation of small protein packing 
i) < 1, the EV approach and the standard Langmuir 
model, eq. (|10[) . are formally equivalent and are allowed 
to be compared, if N = V/(2B 2 ) in the Langmuir pic- 
ture and £' = vq/(2B 2 ). The parameter N, the number 
of fixed binding sites, can then be interpreted as the max- 
imum number of binding spots available for hard spheres 
simply due to packing in the available volume V. In both 
models, Langmuir and EV, the binding constant is refer- 
enced with respect to a standard volume vq = 1/mol. 



E. Numerical evaluation including volume change 

In the Langmuir approach the derivative Q'(x) — 
dQ(x)/dx of eq. (TTTI) is fitted to the experimental data by 
scanning through AH, Kq, and N values until the least 
square deviation (LSD) to the experimental data is min- 
imized. The radius R p is fixed in the Langmuir fittings 
to 1.8 nm. In the EV approach, eq. (fTSj) is solved numer- 
ically and Q'(x), as more generally defined in eq. (|9]) is 
fitted to the experimental data by minimizing the LSD. 
Here the fitting parameters are AH, Kq, and R p . 

One challenge arises because also the gel volume V = 
V(x) is in general a function of x, the protein concentra- 
tion. The electrostatic transfer energy defined by eq. (J3|) , 
in turn, depends on the gel volume. Thus, for separating 
out the electrostatic contributions to binding, the fitting 
of binding isotherms needs take into account the volume 
change. Since a predictive theory for V(x) is out of scope 
of this paper, we satisfy ourselves in the following by 
employing experimental DLS data for R(c s ) and R(x) 
at 7 mM salt concentration as shown in Fig. 2. We fit 
V(x) — 4:n[R(x) 3 — i? 3 oro ]/3 by the empirical function 



(18) R ( X ) = ?(Rn 



Rmin ) 



1 — tanh 



X — Xq 



A 



#min,(20) 



with i? max and i? m j n being the maximum and minimum 
gel radius at x — and x — > oo, respectively, xq is 
the location of the inflection point, and A the distri- 
bution width. For T = 298 K and 7 mM salt we find 
i? max = 172 nm, R min = 129 nm, A = 24000, and x 
can be identified with N(7mM). For the other salt con- 
centrations, where no DLS data are available, V(x) is 
obtained by using the known R ma x = 165.7 nm (17 mM) 
and 157.1 nm (32 mM) (cf. Fig. 2a) and x = N (cf. 
Tables II- IV), while R m [ n = 129 nm is used from the 
7 mM fit. The now only free variable A is employed as 
an additional parameter obtained by least square fitting 
the ITC data. We find A = 31000 (c s = 17 mM) and 
42000 (cs = 32 mM) which appears correlated with the 
change of the sharpness of the binding isotherms N bp (x) 
with c,. 



IV. RESULTS AND DISCUSSION 

A. Gel shrinking by salt and proteins 

In agreement with previous observations) 10 ' 19 ' 25 DLS 
shows that the charged PNIPAM microgel shrinks upon 
the systematic addition of salt as summarized in Fig. 2 
(a) up to a concentration of c s = 1 M. The gel ra- 
dius decreases monotonically and saturates at molar con- 
centrations at a radius close to R ~ 139 nm. The 
best fit of our mechanical balance approach eq. (j6)), also 
shown in Fig. 2 (a), yields very good agreement in al- 
most the whole range of c s . From the fit we obtain 
R(c s — > oo) = Rq = 138.5 nm and a bulk modulus 
K m (Vo) — 204 kPa. Assuming Poisson's ratio to be 1/3, 
valid for neutral polyacrylamide or poly-NiPAm gels ; 71 ' 72 
then Young's modulus E ~ K . Our determined value for 
K is thus fully consistent with recently measured Young 
moduli investigated for varying BIS cross linker density 
at similar temperatures^ In the latter E was measured 
between ~ 80 kPa and ~ 500 kPa for BIS contents of 
2% and 10%, respectively, compared to 5% in this work. 
The good agreement supports our considerations leading 
to eq. ©. 



In Fig. 2(b) we show the measurement of the hydro- 
gel radius R(x) vs. the molar ratio x for the system at 
7 mM salt concentration. Analogously to the salt-only 
case, shrinking of the gel is now observed with increasing 
protein concentration, i.e., with increasing protein load. 
At the highest investigated ratio x ~ 1.1 • 10 5 the radius 
is notably smaller than Rq pointing to binding-related 
network tightening. By fitting to the the empirical func- 
tion (j2"0")l . we find i? m i n = 129 nm. Thus the gel volume 
under full protein load is Vo/V m i n ~ 1.3 times smaller 
than in the neutral reference state. If we assume that the 
charged proteins under full load lead to complete charge 
neutralization and thus TTi on ~ 0, then the bulk modu- 
lus induced by protein binding scales oc (Vo/Vmin) 9 / 4 and 
is ~ 1.7 times larger than without proteins, i.e., the gel 
is roughly 2 times stiffer due to protein sorption. This 
has to be considered a lower bound as it is likely that 
ionic contributions to the osmotic pressure still play a 
role, and maybe protein osmotic effects due to EV inter- 
actions need to be considered. 

However, in order to check to what amount purely elec- 
trostatic effects by proteins induce gel shrinking, we plot 
the description by eqs. ((SJ) and (O also in Fig. 2 (b) 
using the experimental ITC binding isotherm N{,p(x) as 
input. This description is now a prediction and no fitting 
is involved. While the overall shrinking of the gel is rea- 
sonably captured, the model yields a too fast decrease for 
small x. This may point to shortcomings of the mean- 
field cell model (Appendix A), where no ion and protein 
fluctuations are included. For large x, the experimen- 
tal saturation of R m in — 129 nm is not reached indicat- 
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FIG. 2: (a) Microgel radius R{c a ) vs. bulk salt concentration 
c 3 (no proteins). Experimental DLS data (filled circles) are 
fitted by eq. (6) finding Ro = 138.5 nm and a bulk modulus 
K ~ 204 kPa. (b) Microgel radius R(x) vs. lysozyme mo- 
lar ratio x — c l ° % /c m . The green solid line is the theoretical 
prediction from eqs. (6) and (7) based on electrostatic consid- 
erations only. The dashed line depicts a simple empirical fit 
to R(x) = i(i? max - Rmin) [1 - tanh^a] + R mm , [eq. (20)] 
used as input for the fitting of ITC binding isotherms. 



ing that nonelectrostatic effects to gel elasticity play a 
role. Similar unsatisfying performances of simple Don- 
nan models have been observed also in a recent worki^ 5 . 
However, from our comparison it is quite reasonable to 
judge that the dominant effect to gel shrinking by protein 
uptake is of ion osmotic origin. 



B. Characterizing experimental binding isotherms 

1. Langmuir models 

The evaluation of the ITC data by the Langmuir-type 
binding models is presented in Fig. 3 (a). First we com- 
pare the results of different model assumptions for 7 mM 



ionic strength: i) standard Langmuir model vs. ii) ex- 
tended Langmuir (including electrostatic cooperativity) 
with constant volume vs. iii) extended Langmuir includ- 
ing the DLS-measured volume change. From looking at 
the fits by eye and judging from the overall least square 
deviation (LSD) to the ITC data, cf. Tab. II, all fits 
look comparably well and can only be distinguished in 
the presaturation region around x ~ 5 • 10 4 . As discussed 
above fitting is most sensitive to this transition region. 
The fitting parameters are summarized in Tab. II: while 
the heat of binding AH ~ 60 kJ/mol and total number 
of binding sites N ~ 63000 are relatively insensitive to 
the model assumptions, the changes in the binding con- 
stant are big. From the standard model AG = AGo = 
— /3 \sx{K/vq) = —36.6 kJ/mol. In the extended model 
without volume change ii), AGo = —24.5 kJ/mol., i.e., 
more than 12 kJ/mol correspond to ^-dependent electro- 
static contributions. Including the volume change, how- 
ever, has a considerable effect on the electrostatic con- 
tribution which grows by 18 kJ/mol. The latter trend 
is understood by the fact that the monomer charge den- 
sity c g = Ng/V increases with shrinking and the con- 
tributions in AG e ; as given by eq. ([3} rise. Thus con- 
sidering volume changes in charged systems is important 
for quantitative fitting, especially in those systems where 
deswelling is significant. The value of the intrinsic bind- 
ing free energy is AGo = —18.3 kJ/mol for 7 mM ionic 
strength. 

In the next step the extended Langmuir model includ- 
ing volume change has been applied to the other ionic 
strengths 17 and 32 mM as also shown in Fig. 3(a) and 
summarized in Tab. II. Note again that the starting gel 
volumes at x — decrease for increasing ionic strengths, 
see the values for R m ax in Tab. II. Here we first observe 
that the number of Langmuir binding sites N decreases 
with higher ionic strength. The reason is a priori un- 
clear as the Langmuir model assumes a fixed number of 
binding sites independent of ionic strength. We further 
notice that the heat of adsorption AH slightly increases 
with ionic strength. More importantly, however, AGo 
stays fairly independent of c s . Thus, in contrast to the 
standard Langmuir model, the nonspecific electrostatic 
contributions have been consistently separated out and 
the remaining AGo becomes salt concentration indepen- 
dent. On average we find AGo — —18 kJ/mol that might 
be attributable to hydrophobic interactions or possibly 
other local binding effects. 



2. Excluded Volume model 

The fitting of the ITC data by the EV model is pre- 
sented in Fig. 3 (b). Corresponding fitting parameters 
AH, R p , and AGo are summarized in Tab. III. As in 
the Langmuir models the binding isotherms are described 
very well by the fits judging from the overall LSD to the 
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FIG. 3: Fitting of the experimental ITC data (symbols) at 
7, 17, 32 mM salt concentrations by different binding models, 
(a) Langmuir fitting. For 7 mM ionic strength we compare 
the standard Langmuir description (black line), to the Lang- 
muir model including electrostatic cooperativity with con- 
stant gel volume (magenta) and changing gel volume accord- 
ing to eq. (f20|) (red). The data for 17 (green) and 32 mM 
(blue) are fitted using the extended model including the vol- 
ume change. See Table II for obtained fitting parameters, (b) 
Application of the Excluded Volume (EV) binding model to 
the same ITC data as in (a). See Table III for obtained fitting 
parameters. 



experimental data. The results for the heat AH are also 
consistent with the Langmuir fits. This is not surpris- 
ing as this value is determined by the plateau in Q'(x) 
for small x far away from the saturation regime. The 
fitting parameter N is now replaced in the EV model 
by the parameter R p , the effective hard-core radius of 
the protein. The analysis yields values of R p between 
1.8 to 2.0 nm, see Tab. Ill, weakly depending on ionic 
strength. Those numbers are indeed very close to the hy- 
drodynamic radius of lysozyme of about 1.7 irniiS This 
good agreement is actually remarkable and justifies the 
assumptions leading to the EV model, i.e. corroborates 
well with a packing picture of globular proteins. 
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A Go 
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H 
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[kJ/mol] 




7(0 




60100 


61 


2.6-10 b 


-36.6 
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65500 
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1.9-10 4 


-24.5 


89 
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63000 
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1.6-10 3 


-18.3 


44 


17 


165.7 


57500 


62 


2.010 3 


-18.9 


57 


32 


157.1 


42500 


72 


1.310 3 


-17.8 
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TABLE II: Results of fitting to Langmuir models at T — 
298 K. In the 1st row (i) the results of the standard Langmuir 
fit without any correction is shown. In the 2nd row (ii) the 
standard model was corrected for electrostatic cooperativity 
with a constant gel volume. In rows 3 to 5 the change of gel 
volume was considered in the fit according to eq. (|20|) . LSD 
denotes the least square deviation to the experimental data. 
Lower values correspond to better quality of the fits. 
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The magnitude of the intrinsic adsorption energy 
|A£ro| in Tab. Ill are between 15 and 18 kJ/mol match- 
ing closely the ones from the Langmuir model in Tab. II. 
A very small salt dependence of AGo remains indicat- 
ing a slightly less accurate subtraction of the nonspe- 
cific effects in this model. However, the salt concen- 
tration dependence is quite small and on average we 
find AGo — — 17 kJ/mol that has to be attributed to 
specific local binding effects. In Tab. Ill we also show 
the results of EV model fitting to experimental data 
gathered at a higher temperature, T — 303 K. The in- 
trinsic binding affinity increases slightly to an average 
AGo — -18 kJ/mol which corroborate with the typi- 
cal thermodynamic signature of increasing hydrophobic 
association at enhanced temperatures^ 
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1.79 


70 
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65 


-19.6 
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TABLE III: Results of fitting to the Excluded Volume (EV) 
model at T = 298 K (top) and T — 303 K (bottom). The 
change of gel volume was considered in the fit according to 
eq. (]20p . LSD denotes the least square deviation to the ex- 
perimental data. Lower values correspond to better quality 
of the fits. 



3. Donnan potentials A(j>(x) and the binding constant K(x) 

The Donnan potential A(j>(x), the electrostatic binding 
energy AG e i(x), and so the total binding constant K{x) 
depend on protein load, and thus, in turn, on the molar 
ratio x. The Donnan potential A<p(x) as resulting from 
the fitting to the extended Langmuir model is plotted in 




5.0xl0 4 LOxlO 5 1.5xl0 5 



x 

FIG. 4: (a) Mean (Donnan) potential A(j>(x) inside the gel 
vs. molar ratio x according to the Donnan form eq. (1) and 
(4) obtained from fitting of the extended Langmuir model to 
the ITC data. For 7 mM we compare the results from fitting 
with constant volume vs. variable volume, (b) Total binding 
affinities — In K(x)/vo (solid lines) and total electrostatic con- 
tribution AG e ;(a;) (dashed lines), cf. eq. (3), to the protein 
transfer energy from the extended Langmuir model for the 
different ionic strengths 7, 17, 32 mM, as labeled. Also plot- 
ted is the hard-sphere Carnahan-Starling contribution fics 
from the EV model fit at 7 mM (dashed-dotted line). 



Fig. 4(a) for all ionic strengths. In the case of c s = 7 mM, 
at the beginning of the titration, x ~ 0, the potential 
magnitude is a considerable 1.77 fc^T/e (~ 44 mV). 
The potential quickly decreases with protein load due 
to charge neutralization, cf. eq. and saturates for 
x > N due to saturation of Nj,p(x). The decrease of gel 
volume in the adsorption process has a notable effect on 
the potential at large x due to the accompanying increase 
in monomer charge density. Increasing the ionic strength, 
on the other hand, considerably lowers the potential at 
small load, while surprisingly A<f> remains roughly inde- 
pendent of c s for high load. The reason for the latter 
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is that less proteins bind for higher ionic strengths, i.e., 
there is less charge neutralization. The Donnan potential 
is expected to be highly correlated with the surface po- 
tentials of the microgels which govern colloidal stability 
in solution^ 

The total electrostatic contribution AG e i(x) to the 
binding constant is plotted in Fig. 4(b) together with the 
total binding free energy ~kgT\ii[K(x)/v ]. The curves 
AG e i(i) and — fcsTln[A(x)/i>o] are parallel with an off- 
set given by the intrinsic part AGo , as defined in eq. (H"3l) . 
The electrostatic contribution for very small load (x = 0) 
is big and contributing a favorable 15 fc^T (~ 37 kJ/mol) 
for 7 mM and 6 k B T (~ 15 kJ/mol) for 32 mM to the 
total binding affinity of about 22 fcgT and 14 fcgT, re- 
spectively, at the beginning of the titration. We find that 
the Born contribution (second term in ([3])) constitutes a 
significant (favorable) 2-3 ksT to the total electrostatic 
energy and is thus not negligible. All curves decrease 
in their absolute values and saturate for values x > N, 
when the molar ratio exceeds the number of binding sites 
in the Langmuir model. Note that for 7 mM the satura- 
tion value of K (x) is close to the value of K obtained from 
standard Langmuir fitting, c.f., Tab. II, because fitting is 
most sensitive to the presaturation region x ~ A above 
which K(x) remains constant. This statement holds also 
for the other salt concentrations when K{x) is compared 
to previous standard Langmuir fitting^ 

Also shown in Fig. 4 (b) is the Carnahan-Starling con- 
tribution fics m the EV model. This entropic penalty 
due to hard sphere packing slowly rises with x in the 
regime x < A until quickly increasing to unfavorable 3 
ksT in the saturation regime x > A. 



C. Interpretation of Langmuir vs. EV model 
results 

In section lTl! D 4l we argued that the EV model is equiv- 
alent to the Langmuir approach in the low packing regime 
(rj <C 1), if A, the number of binding sites, is equated to 
V I (2B2) the number of free spots just by packing. Indeed 
if one now considers a mean gel radius of R — 150 nm 
and protein diameter a p = 3.8 nm obtained from fit- 
ting above, we end up with A = V/(2B 2 ) ^ 59000, to- 
tally consistent with numbers obtained from fitting to 
the Langmuir models above, cf. Table II. This agree- 
ment implies that the strict Langmuir assumption of a 
fixed set of binding sites can still be considered an inter- 
pretable quantity, even if the nature of the bound state is 
ill-defined. From that point of view, the decrease in A in 
the Langmuir models with increasing ionic strengths can 
be understood: for increasing c s , the gel volume V de- 
creases and packing penalties are becoming more impor- 
tant for a smaller number Nb p of bound protein. Hence, 
in standard Langmuir fitting of ITC data for sorption 
to soft materials,— ^ 20 ' 22 ' 23 ' 31 where it is likely that not 
a condensation-like binding of proteins takes place, the 
stochiometry of binding may be interpretable by packing 



effects. 

As we demonstrated in Fig. 4 the electrostatic energies 
and therefore the binding affinity K{x) are functions of 
the molar ratio x. How do we interpret a constant K as 
deduced from standard Langmuir fitting?— i L 20 i 22 i 23 i 31 As 
argued in the Methods section fitting is most sensitive 
to the presaturation region x ~ A. Indeed the data in 
Fig. 4 (b) shows that K{x ~ A) equals the values of stan- 
dard Langmuir fitting, see Tab. II. Thus, we can conclude 
that a binding constant obtained from a standard Lang- 
muir fit is a reasonable number which can be interpreted 
as binding affinities in the presaturation regime, i.e., in 
the intermediate to high protein load regime, where also 
volume changes are not so large anymore. However, for 
small x <C A (small protein load) our separation into 
electrostatic and hydrophobic contributions shows that 
binding affinities can be about 1000 times larger than in 
the presaturation regime. This may have implications 
for the modeling and interpretation of protein binding 
kineticsiiSiS^ 

It may be on first glance remarkable that in both 
models, extended Langmuir and EV, the consistent sep- 
aration of electrostatic and hydrophobic effects yields 
the same number for the hydrophobic (intrinsic) bind- 
ing affinity. In our system we find AGo — —7 k B T which 
must be attributed to nonclcctrostatic binding effects and 
constitutes roughly 1/3 or 1/2 of the total binding affin- 
ity in the small and high load regimes, respectively. The 
very weak dependence of AGo with salt concentration in 
both models indicates a successful separation of nonspe- 
cific electrostatic and intrinsic effects in our treatment. 
However, note that the agreement can only be established 
if the nature of the bound state is identically defined in 
both models as discussed in section IIIID4I Since at low 
packing we exemplified total equivalence of both models, 
the agreement for AGo is thus not unexpected. Prob- 
ably for higher packing densities than observed in this 
study (where r/ < 0.18), the results from both models 
may disagree stronger. Thus for not too high packing, 
one model does not seem superior over the other in de- 
scribing experimentally measured binding isotherms, if 
correctly interpreted. 

Furthermore we would like to comment on the magni- 
tude of AGo — ~7 ksT for the instrinsic interaction 
of lysozyme with the PNIPAM network. If methane- 
methane interactions are taken as reference with attrac- 
tions on the order of 2-3 fc^T, then 7 ksT correspond 
to 2-3 hydrophobic protein-PNIPAM contacts on aver- 
age which seems reasonable. Increasing temperature led 
to an increase in AGo conform with the signature of hy- 
drophobic interactions^ A recent study on a similar sys- 
tem showed hardly uptake of lysozyme by an uncharged 
PNIPAM microgelj 2 ^ Reasons for this discrepancy may 
be the different gels preparation state, e.g., larger pore 
sizes within the gel. Alternatively, maybe we overesti- 
mate the effects of hydrophobicity and other local effects, 
such as salt bridges, play a more important role as ex- 
pected. 
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Finally, we note that effective net charge of chicken 
egg white lysozyme as used in this work may be slightly 
larger on average due to the protonation effects (z p ~ 7 
to 8) within the gelj^ 3 - However, using z p = +7.5 or +8 
in our analysis we end up with a similar A Go ~ — 7 fc^T; 
the reason is that while the prefactor in the electrostatic 
contribution (3) rises, the Donnan potential (2) decreases 
quicker with load. These effects roughly cancel each 
other for our particular system. 



arise, however, if the systems of practical interest exhibit 
microscopic irreversibiltiy and the final protein binding 
cannot be considered in equilibrium terms^^ 
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V. CONCLUDING REMARKS 

In summary we have systematically introduced bind- 
ing models to characterize the physical interactions in 
the process of equilibrium protein sorption to microgels. 
The binding models separate out electrostatic coopera- 
tivity and include deswelling effects. The analysis yields 
a hydrophobic, salt-independent interaction of around 
A Go ~ — 7 ksT for our core-shell system and hen-egg 
white lysozyme. 

The found A Go constitutes roughly 1/3 or 1/2 of 
the total binding affinity in the small and high load 
regimes, respectively. For small protein load, binding 
affinities can be about 1000 times larger than in the pre- 
saturation regime due to high net charge of the protein 
which is gradually neutralized upon binding. This may 
have implications for the modeling of protein binding at 
early stages of sorption, especially to understand kinetic 
rates i 19 ! 75 " — Moreover, the Donnan potential is expected 
to be highly correlated with the surface potentials of the 
microgels which govern colloidal microgel stability in so- 
lution^ 

We find gel deswelling to be mostly of electrostatic ori- 
gin. However, the gel becomes at least 2 times stiffer at 
high load pointing either to more specific effects at high 
load (e.g., cross links) or electrostatic correlations not ac- 
counted for in mean-field PB cell model approaches. The 
change and control of material properties upon protein 
load is essential for functionality, — and suggest challeng- 
ing investigations in future. 

In two complementary models, extended Langmuir and 
EV, the consistent separation of electrostatic and hy- 
drophobic effects yields the same number for the hy- 
drophobic (intrinsic) binding affinity. The agreement can 
only be established if the nature of the bound state is 
identically defined in both models. Details of the bound 
state can be inferred, for instance, from small angle scat- 
tering^ We work at low protein packing; hence, one 
model does not seem superior over the other in describing 
experimentally measured binding isotherms, at least for 
not too high loads (packing). With that a more quan- 
titative interpretation of binding data iLLSLSSiS&SiiiSSjSL 
in terms of separate physical interactions is possible in 
future. Especially we have shown that fitting based on 
standard Langmuir models yields intcrprctablc binding 
affinities and stochiometryi£i 7 - i 20 ' 22 i 23 ' 31 More challenges 



Appendix A: Linearized Poisson-Boltzmann (LPB) 
cell model 

To obtain the electrostatic contributions to the trans- 
fer of charged hard spheres from bulk to gel we make use 
of the Poisson-Boltzmann cell model. While often nu- 
merical solutions have been employed ^ r 38 1 40-43,51,78 f or 
weak perturbations the linearized form can be treated 
analytically! 33 ' 39 

In the cell model each protein is represented by a ho- 
mogeneously charged sphere with radius R p and valency 
z p centered in a spherical cell with radius R c and vol- 
ume V c . The gel is assumed to be made up of Nb p such 
cells whose dimensions are determined by the number of 
proteins in the volume, i.e., 

where Nb p is the total number of proteins in the vol- 
ume V. Each cell is in contact with a reservoir with 
monovalent salt concentration c s . Inside the gel each cell 
additionally contains a fixed number of iV c = N g /Nb p 
charges from the charged network monomers. The 
mean monomer charge concentration in one cell is thus 
N c /V c = N g /V — c g . Each cell must be electroneutral 
on average, i.e., in the gel it holds 

c s e- A * - c s e A * + z gCg + z p N bp /V = 0, (A2) 

where we neglected the vanishingly small protein con- 
centration outside of the gel. For high protein load, R c 
becomes comparable to R p and the cell volume needs in 
principle to be corrected by the protein volume. However, 
for our systems at highest protein load (R p /R c ) 3 < 0.1, 
and the correction is negligible for small and intermedi- 
ate loads. The solution of (IA2[) is the modified Donnan 
potential 

A0 = bx[v + VlF+l] (A3) 

with y = izgCg + z p Nb p /V)/(2c s ) which describes the 
difference in the mean electrostatic potential with respect 
to the bulk reference state where we set <\> = 0. In the 
bulk the protein concentration is typically vanishingly 
small, thus electroneutrality dictates c + = c_ = c s to a 
very good approximation. 
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We now focus on the proteins inside the gel. Here 
we assume the PE network in the gel to be flexible and 
fluid-like and thus the N c charged monomers to behave 
like mobile counterions to the protein. The PB equation 
in spherical coordinates is then 



~W0" 

r 



-4irl B [(z g cexp(-z g 4>i) (A4) 
c s exp(A</> + 4>x) + c s cxp(-A^ - cf>i)], 



where we made the ansatz <p(r) = A(f> + <pi(r), with A<p 
being the constant mean (modified Donnan) potential in 
eq. (|A3|) and 4>\ being a perturbation induced by the 
protein. Note that if the monomers were assumed to 
be just a fixed, homogeneous background, they would 
not couple to the field and cexp(— z g <pi) needed to be 
replaced by the fixed concentration c g . The constant c 
is defined by conservation of the number of monomer 
charges in the cell via 



N c = c d J rexp(-z s ^ 1 ). 



(A5) 



It also holds that the average potential equals the Donnan 
potential, i.e, 



A</> = — / d 3 rif>(r 



(A6) 



In the limit of large cell sizes, i.e., R c — > oo (or 
Nb p — > 0), it follows that quickly B — > 0, and A — > 
z p l b exp(Ati? p )/(l + kR p ). The LPB solution simplifies to 



<j>{r) = A(j>- z p N bp /(Vc in ) + - 



z p Ib 



,(A9) 



where we kept the important leading order term in 1/V. 
An illustrative sketch of this potential distribution in the 
cell model is given in Fig. [5j A(j> is the average poten- 
tial, while for small protein load the potential at the cell 
boundary 4>(R C ) — A(f> — z p N bp /(Vci n ) and at the parti- 
cle surface (j>(R p ) ~ A<j>- z p N bp /(Vc m ) + (z p l B )/[R P (l + 
kR p )]. For not too high salt concentrations c s < c g , it fol- 
lows that the Donnan potential A<p > 1, and the internal 
salt concentration can be well represented by q„ ~ 2c g , 
as the coion concentration (oc exp[— |A0]|) in the gel be- 
comes negligibly small. In the bulk solution in the dilute 
protein limit analogously to eq. (|A9|) we find 



0(r) 



zJi 



-K b (r-R. p ) 



(A10) 



1 + K b Rp r 
where n b = ^/8ttIbc s is the bulk inverse screening length. 



For not too large protein charges and thus (j>i <§; 1, we 
can linearize the exponentials in the PB eq. (|A4|) with 
respect to <j>i which yields 

-(r4>)" = 4nl B z p N bp /V + k 2 (4> - A$), (A7) 
r 

where we have identified the internal (gel) charge den- 
sity — z p N bp /V = z g c g + c s exp(— A<fi) — c s exp(A</)) from 
conditions (A5) and (A6), the internal concentration 
Cin = Cg + c s exp(— Atfi) + c s exp(A</)), and the internal 
inverse screening length n = y/inlBCm- The LPB equa- 
tion can be solved with respect to the boundary condi- 
tion that the electrical field on the sphere surface is fixed 



by cf>'(R p ) = 
ary $'(R C ) 
work o 33 ' 39 ) 



Zjyl 



i/R p and it vanishes at the cell bound- 



0. The final solution is (see also previous 



(r) = Aif> + Mr) 



c~ c 

A0 - z p N bp /{Vc in ) + A + B — 

r r 



with constants 



A = 
and 

B = 



z p l B e KR " r 
1 + kR v 



1 _ e -2K(R a -R p ) 



{ K R p - 1) (kRc - 1) 



-i -l 



(kR p + 1) (kR c + 1) 



Zn I 



plB 



1 + K,Rp 



, K (2R C -R p ) ( kR c - 1) _ eKRp {nRp - 1) 

(kR c + 1) 



(kR p + 1) 



A((> 




Y R 




FIG. 5: Illustrative ID sketch of the electrostatic potential 
distribution in the cell model: the charged proteins (blue 

(A8) 

spheres) with radius R p sit equidistantly in the gel in their 
cells with radius R c . 4>(R P ) is the potential at the protein sur- 
face, <j>(Rc) the potential at the cell surface, and A(f> the mean 
(Donnan) potential. In our approach <f>i(r) is a perturbation 
of <j)(r) around A<f>. 



The transfer (Gibbs) free energy (or chemical poten- 
tial) to bring a protein from bulk solution into the gel 
can then be calculated by the difference of work of 
charging the sphere against the surface potential (j)(R p ), 
/3AG = J 2p &z p [4> g {R p ) - 4> b (R p )) in a cell in the gel (g) 
vs. bulk (b). We obtain the leading order contributions 
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(in the low salt and low protein limit) 

zlN bp 



/3AG, 



ZpAtfi - 



ZplB 



2 -R E 



2c g V 

KgRp 
+ KgRp 



(All) 



1 + KbRp 



where we defined K g = ^8irl B Cg. We did not explicitly 
integrate over z p in the second term in (|A9[) as it is a con- 
stant background contribution not immediately involved 
in the charging process of one particle. It is notewor- 
thy that we obtain the same functional form for AG e ; if 
the monomer charges are not assumed to be mobile, al- 
beit with a \/2 smaller internal inverse screening length 
yj AirlBCg and a factor of 2 in front of the second term. 

By a one-to-one comparison of the leading order ex- 
pression (|A11I) to the result from employing the full ex- 
pression (|A8I) we find that the error in AG e i is less than 
fcflT over the whole range of molar ratios and salt concen- 
trations used in this work. By detailed inspection of the 
behavior of (|A9I) we observe that the reason of the sur- 
prising accuracy is a fortuitous cancellation of errors of 
higher order terms at large N bp . This fact may shed some 
doubt on the general applicability of simplified eq. (jAllj) 
but note that our parameters (protein valency, salt con- 
centrations, monomer charge densities, etc.) are typical 
for a wide variety of experimental systems. However, 
in general the PB approach is expected to break down 
for very high protein valencies and small proteins, when 
\<pi\ ^> 1, and strong Coulomb correlations play a role.— 

The free energy (|A1 1|) above considers only ionic con- 
tributions to solvation of a fixed lattice of spheres, i.e., 
it neglects the electrostatic contributions from the in- 
teraction between the proteins, i.e., the energy penalty 
of overlapping double layers. However, in the fluid-like 
hydrogel protein matrix it is reasonable to assume that 
proteins can wiggle or move around and are not rigidly 
fixed to lattice positions. Due to such fluctuations the av- 
erage surface potential will actually be higher than given 
in (|A8I) . To estimate the interaction contribution we look 
at the expansion of the excess chemical potential in terms 
of virial coefficients, i.c, in first order 



where Bo 



0H = 2B 2 /V C , (A12) 
\ J d 3 r [exp(-0W(r)] - 1) and 



W(r) = Whs 



zIIb 



D -n(r-R p ) 



1 + kR p r 



(A13) 



is the protein-protein interaction potential, split up into 
hard-sphere and the Debye-Hiickel contribution accord- 
ing to (IA8|) . Using condition (| A6|) . linearizing the ex- 
ponent in the defining equation for B2, and splitting the 
chemical potential into hard-sphere and electrostatic con- 
tributions we find 



0H = Phhs + 



ZpN bp 
2^V 



(A14) 



Thus in leading order the electrostatic protein-protein in- 
teraction contribution exactly cancels the second term in 
(lAlip and we end up with the final result for the electro- 
static transfer free energy 



fiAGel = ZpA(j) - 



Zpls 

2R n 



KgRp 



1 



KgRp 



KbRp 
1 + KbRp ) 



(A15) 



It is interesting to note on the simplicity of eq. (|A15|) 
which describes naively the transfer of charge z p into the 
average potential A</> and, in the second term, the dif- 
ference in Born solvation free energies in a homogeneous 
medium with salt concentrations c s and c g . Thus, the in- 
homogeneities introduced by the cell model assumption 
as depicted in Fig. 5 cancel out (in linearized theory) 
if particle fluctuations are allowed, and the naive form 
(|Al5l) holds. 

The estimation of ionic osmotic contribution to pres- 
sure in the presence of proteins seems less simple. Micro- 
configuration of proteins induce inhomogeneities and the 
pressure is not anymore given by the simple expression 
for 7r in eq. ([2J . The ion contribution to the osmotic pres- 
sure is usually estimated from the ionic concentration at 
the cell surface (in the cell model) where the electrostatic 
pressure vanishes r^r— i2&2iLZ& 

7rf on ~ c s exp(^(i? c )) + c s exp(-4>(R c )) - 2c s . (A16) 

Expanding to 1st order in </>i, we obtain 

7rf on ~ 7T ion + 2c s 0i {R c ) sinh(A^), (A17) 

which recovers 7rf on ~ ~Ki on in the limit for vanishing pro- 
tein concentrations and A0(y) is the protein-corrected 
Donnan potential (|A3I) withy = (z g c g + z p Nbp/V)/(2c s ). 
Note that this expression does not consider fluctuations 
in protein positions which is likely to be an important 
effect to consider in future studies. 



Appendix B: The standard Langmuir model in the 
canonical ensemble 



Consider a finite region in space with N identical and 
independent binding sites available. We denote the num- 
ber of bound proteins by Nb p and define the fraction of 
bound particles by = Nb p /N. The number of available 
binding states is the n 29 ' 30 



W 



N b p\(N-N bp )V 



(Bl) 



from the combinatorial possibilities of distributing Nb p 
indistinguishable particles on N sites, and £ is the par- 
tition sum of a single particle in the bound state. The 
Boltzmann entropy is defined by 



^- = hxW 

k B 



(B2) 
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leading (within a constant) to the entropy per binding 
site 

= -e in e - (i - e) in(i - e) + e m( va 3 XB3) 

where we defined C t Nbp — (vo/A 3 ) Nbp in terms of an effec- 
tive configurational volume vq divided by the cubed ther- 
mal (de Broglie) wavelength A 3 . 'Effective' means that 
also restrictions on vibrational and orientational degrees 
of freedom are adsorbed in the number vo, not purely 
translational effects if it all. The (canonical) Hclmholtz 
free energy of the system is 

f3F = (3Fid - 4- + (3N bp AG, (B4) 
Kb 

where we introduced the canonical ideal gas free energy 
/3F id = (n- N bp )[ln((n - N bp )A 3 /V) - 1] with total par- 
ticle number n. Then (n — N bp )/V is the density of un- 
bound particles in a total volume V (assumed to be much 
larger than the binding region), and the adsorption free 
energy AG associated with the binding of one protein 
to one Langmuir site. The free energy / = (3F/N per 



binding site is then 

/ = jf(n~ N bp )[ln((n - N bp )A 3 /V) - 1] (B5) 
+ ein6 + (l-6)ln(l-e) 
- ein(VA 3 )+/39AG. 

The minimization of the free energy with respect to 
the number of bound protein df /dN bp = yields then 
the final relations for the fraction of bound particles 
in dependence of the (unbound) particle concentration 
c p = (n— N bp )/V. We obtain the final result 

K = cxp(-/3 AG)v = 9 . (B6) 
(1 - 6)c p 

The standard volume t>o depends on the exact nature of 
the bound state and is typically not known. Thus the 
prediction of absolute binding free energies is difficult. 
In literature often the standard volume l/mol~ 1.6 nm 3 
is employed. This is not unreasonable if one assumes that 
still rotational and vibrational modes in the bound state 
take place on molecular, i.e, nanometer scales. 
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